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, Abstract 

A variable-skew oblique wing offers a sub- 
stantial aerodynamic performance advantage for 
aircraft missions that require both high effi- 
ciency in subsonic flight and supersonic dash or 
cruise. The most obvious characteristic of the 
oblique-wing concept is the asymmetry associated 
with wing-skew angle which results in significant 
aerodynamic and inertial cross-coupling between 
the aircraft longitudinal and lateral-directional 
axes. This paper presents a technique for synthe- 
sizing a decoupling controller while providing the 
desired stability augmentation. 

The proposed synthesis procedure uses the con- 
cept of explicit model following. Linear quad- 
ratic optimization techniques are used to design 
the linear feedback system. The effectiveness of 
the control laws developed in achieving the desired 
decoupling is illustrated for a given flight con- 
dition by application to linearized equations of 
motion, and also to the nonlinear equations of 
six degrees of freedom of motion with nonlinear 
aerodynamic data. 

Nomencl ature 

A,B system matrices, 

h altitude 

J cost functional 

K gain matrix 

M Mach number 

OWRA oblique-wing research aircraft 

p roll rate , deg/sec 

Q,D state weighting matrix 

q pitch rate, deg/sec 

R control weighting matrix 

r yaw rate , deg/sec 

s complex frequency 

t time 


♦Associate Professor, Electrical Engineering 
Dept., Member AIAA. 

♦♦Aerospace Engineer. 

^Aerospace Engineer, Member AIAA. 


Cal ifornia 

u input vector 

v velocity, ft/sec 

x state vector 

a angle of attack, deg 

6 sidesl ip angle, deg 

6 control surface deflection 

6 pitch angle, deg 

A wing skew angle, deg 

<)> bank angle, deg 

Subscripts 

aL left aileron 

aR right aileron 

hL left horizontal tail 

hR right horizontal tail 

m model 

p aircraft 

Superscripts 

m number of inputs 

n number of states 

T matrix transpose 

Introduction 

The advantages of an oblique wing were first 
noted in the 1940's. However, not until recently 
have the interest, technology, and mission of an 
oblique-wing design evolved into a full-scale 
flight research program. Dryden Flight Research 
Facility of NASA Ames Research Center and the 
U.S. Navy are developing an oblique-wing research 
aircraft (OWRA). Gregory,! and Wiler and White2 
have outlined potential advantages and disadvan- 
tages of this type of airplane. Theoretical and 
wind tunnel studies have shown that a variable- 
skew oblique wing offers a substantial aerodynami 
performance advantage for aircraft missions that 
require both high efficiency in subsonic flight 
and supersonic dash or cruise. 


The most obvious disadvantage of the oblique- 
winy concept is the asymmetry associated with 
wing-skew angle. This asymmetry results in sig- 
nificant aerodynamic and inertial cross-coupl ing 
between the aircraft longitudinal and lateral- 
directional axes. 


as the model would respond to a given pilot com- 
mand. More precisely, the model -fol 1 owi ng problem 
can be stated as follows. 

Given the linearized aircraft dynamics, 


The test bed for OWRA will be the NASA F-8 
digital fly-by-wire aircraft. This aircraft will 
be modified by the removal of the current high 
wing. and installation of a composite wing-and- 
pivot assembly. A major part of the OWRA program 
will be the synthesis of a flight control system 
which will provide acceptable stabilization and 
decoupling across the Mach, angle-of-attack , and 
wing-skew envelope. These aircraft stability and 
decoupling requirements are ideally suited for the 
application of modern control theory techniques to 
the solution of problens associated with OWRA. 

The potential advantages of the oblique-wing 
concept can only be realized by development of 
related technologies. Foremost among these is the 
control system architecture to compensate for 
cross-coupled aircraft responses, while presenting 
the crew with the feel of a conventional airplane. 
Current typical design procedures synthesize 
aircraft controllers based on solutions of two, 
or at most, three degrees of freedom. However, 
the added OWRA stabilisation and decoupling prob- 
lems require at least ive degrees of freedom 
simultaneously. 

Model following ha:: been a popular method for 
the design of multi var able control systems, and 
shown to be amenable tn the solution of many air- 
craft control problems.. Commencing with the work 
of Kalman 3 and Tyler ,4 extensive use has been made 
of the linear optimal control theory in design of 
model -fol lowing controllers. Yore 5 used this 
method for simultaneous stability augmentation 
and mode decoupling. While optimal control the- 
ory provides an extremely flexible synthesis tech- 
nique, various structural approaches have been 
adopted in synthesizing model -foil owing control- 
lers. 5 A controller based on perfect model fol- 
lowing (such as perfect matching of the dynamics 
of the compensated plant to those of the model) 

was presented by Alag, Kempel , and Pahle^ for 
control of OWRA. Though the desired degree of 
decoupling was achieved, the control surface 
activity was excessive and alternate controller 
development was required. 


x p = ApXp + BpUp ( 1 ) 

where x p eR n , UpcR 171 , A p , and B p are matrices of 
appropriate dimensions, find the control u p • 
such that the aircraft states Xp approximate 
“reasonably well" model state vector x m . The 
model state vector xm is defined by the equation: 


x m = Am x m + B m Um 


( 2 ) 


where x m eR n , u m eR m , Am, and Bm are matrices of 
appropriate dimensions. 


For OWRA, the state and input vectors are 
given by 


x 



velocity, ft/ sec 
angle-of-attack, deg 
sidesl ip angle, deg 
bank angle, deg 
pitch angle, deg 
roll rate, deg/sec 
pitch rate, deg/sec 
yaw rate, deg/ sec 


u = 


6 hi." 
6 hR 
«aL 
«aR 
_ S R . 


left horizontal tail deflection, deg 
right horizontal tail deflection, deg 
left aileron deflection, deg 
right aileron deflection, deg 
rudder deflection, deg 


The desired model of the ai rcraft , defined by 
matrices ^ and Bm, as well as the aircraft mat- 
rices Ap and Bp, are given in Table 1. The desired 
model used in this study is a modification of the 
zero-wing-skew configuration at the same flight 
condition. The aircraft matrices correspond to a 
flight condition of Mach 0.8 and -an altitude of 
20,000 ft at 45° wing skew. The elements of Am 

were modified to increase the short period and 
dutch-roll mode damping, and to provide improved 
roll and spiral mode characteristics. 


This paper presents the use of models in 
design of linear feedback systems by means of 
linear-quadratic optimization. The linear control 
law developed provides the decoupling as well as 
the desired stability c.ugmentation . The effec- 
tiveness of the control law is illustrated by time 
responses from linearized equations of motion, and 
nonlinear equations of six degrees of freedom of 
aircraft motion for a c.iven flight condition. 

Model -Fol 1 owing Systems 


Explicit Model -Following Systems 

There are two configurations of model fol- 
lowing, known as implicit model following and 
explicit model following. As Fig. 1 shows, in 
implicit model following, the control inputs to 
the plant are formed from the aircraft states 
and pilot input. No dynamic coupling exists 
between the model states and the closed-loop 
plant; the model state x^ appears only in the 
performance index. 


Problem Definition 

The concept of model following is useful when 
an ideal set of plant equations of motion can be 
specified. The objective of model -fol lowing 
flight control is to force the aircraft to respond 


Figure 2 illustrates the explicit model 
following in which the model states must be 
generated for use in forming the control input. 
Explicit model following requires the simulation 
of the model as a part of feedforward controller. 
Alignment of plant and model in the presence of 
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uncertainties, such as unknown parameters and ran- 
dom disturbances, requires this type of control. 
This enables a continuous correction of the errors 
between model and plant states even in the pres- 
ence of unknown disturbances. 

The dynamics of the plant and model are 
governed by the following linear state equations: 

*p = ApXp + BpUp . (1) 

x m = AmXfli + Bm u m (2) 

Defining an augmented state vector x, the dynamics 
of the system are expressed as 

x = Ax + Bup (3) 

where 



The pilot input u m is modeled as a constant, 
an assumption that does nut overly distort 
the reality of the situation and allows a com- 
plete analysis of the problem from a theoret- 
ical viewpoint. 

In explicit model following, tne control func- 
tion up is required to minimize a performance 
index given by 

- DO 

J = j J Q C(x p - Xm) T q(x p - xm) + UpRu p ] dt (5) 

where Q ^ Q, R >0 are weighting matrices. Using 

the augmented state vector x, the index J can be 
rewritten as 


response of pitch and yaw angular rate and bank 
angle to an. elevator command input. Significant 
yaw rate and bank angle are generated as a result 
of pitch command. Of particular interest is the 
large change in bank angle, indicating a high 
degree of cross-coupling. 

Table 2 gives the gain matrices for the 
explicit model -foil owing configuration. Figure 4 
indicates the linear closed-loop system response 
for the command input with the explicit model - 
following gains. The aircraft response is indi- 
cated by the solid line, and the model response 
by the dashed line in Figs 3 and 4. The model- 
following response is considered satisfactory with 
considerable attenuation in the degree of coupling, 
as indicated by the bank angle response of the air- 
craft. The model response in bank angle is zero. 
Figure 5 indicates the closed-loop response for 
the same case but with the aircraft now repre- 
sented by a nonlinear simulation of 6 degrees of 
freedom for the same flight condition. The non- 
linearities include the constraints on the con- 
trol surface rate and position. 

Figure 6 indicates the response of the open- 
loop system to 1° aileron command input. Bank 
angle, yaw angular rate, and pitch angular rate 
are shown. The relative pitch coupling is not 
as severe for this case as the roll coupling is 
for the elevator command; however, coupling is 
still evident. 

Figure 7 indicates the linear closed-loop 
system response to the same aileron command. A 
dashed line denotes the model response, and a 
solid line, the aircraft response. The relatively 
small coupling in the pitch axis was not con- 
sidered severe, as indicated by the aircraft 
pitch-angular rate response. The model pitch- 
angular rate response was zero. The objective 
in this case was to provide adequate lateral - 
directional dynamics. Figure 8 shows the closed- 
loop response for the same case with aircraft non- 
linear simulation. 



If pair (A,B) is stabilizable and pair (A,D with 
D t D = Q) detectable, the optimal control up which 
minimizes J is given by 

Up = Kx = K X pXp + KxmXm + KumUm (8) 

The model equations must be simulated as a part of 
the feedforward controller because the controller 
requires model states. 


Conclusions 

The method presented describes a decoupled 
control for a highly coupled unsymmetric aircraft. 
The method develops an 'explicit model -fol lowing 
control law by means of linear quadratic optimi- 
zation techniques. The results indicate that 
the method does obtain the decoupling incorpo- 
rated in the ideal model for the flight condi- 
tion considered. 

Evaluation of the control system on nonlinear 
equations of six degrees of freedom of motion for 
the flight condition considered shows promising 
results for implementation as a candidate control 
system for the aircraft. Work is in progress to 
investigate gain scheduling requirements and to 
obtain piloted simulation results. 


Results 

The degree of coupling in the open-loop air- 
craft is illustrated by response to application 
of a 1° command input at 1 sec and returned to 
zero at 3 sec (either by elevator or aileron 
input). Figure 3 illustrates the open-loop system 
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TABLE 1. - AIRCRAFT AND MODEL MATRICES* 


(a) Aircraft matrices 



r o.oo9f- 

22.0707 

10.5479 

-0.1341 

-32.1127 

0.0057 

-0.0005 

-0.02651 



-0.0001 

-0.7826 

0.0958 

0.0000 

0.0000 

0.0030 

0.9926 

-0.0003 



0.0000 

-0.0592 

-0.2908 

0.0387 

-0.0002 

0.0259 

0.0001 

-0.9920 



0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0247 


A p = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 



0.0000 

33.1432 

-53.6933 

0.0000 

0.0000 

-3.1250 

2.0552 

1.7210 



-o.ooo;: 

-8.6816 

0.7975 

0.0000 

0.0000 

0.1679 

-1.0352 

0.1810 



L 0.0000 

-1.0092 

10.7521 

0.0000 

0.0000 

-0.0213 

0.0080 

-0.7129J 



f 0.0844 

-0.0309 

-0.2210 

-1.2572 

3.6598^ 






-0.0974 

-0.0974 

-0.0198 

-0.0302 

0.0000 






-0.0166 

0.0166 

0.0008 

-0.0005 

0.0647 






0.0000 

0.0000 

0.0000 

0.0000 

0.0000 





Bp - 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 






12.9804 

-22.2654 

15.8467 

-11.8422 

13.2774 






-9.407:t 

-10.8555 

-1.2311 

0.8797 

0.5694 






L 1.9854 

-2.2579 

0.5262 

-0.3276 

-6.2493J 








(b) 

Model matrices 






r-o.oo7' 

23.5966 

0.0000 

0.0000 

-32.1129 

0.0000 

0.0000 

0.0000 

1 


-o.ooo:. 

-1.1CS2 

0.0000 

0.0000 

0.0000 

0.0000 

0.9909 

0.0000 


0.0000 

0.0000 

-0.2852 

0.9387 

0.0000 

-0.0148 

0.0000 

-0.9919 


0.0000 

0.0000 

.0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

-0.0120 

A m " 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 


0.0000 

OiOQOO 

-44.3777 

0.0000 

0.0000 

-4.5000 

0.0000 

14.6000 



0.0001 

-12.1514 

0.0000 

0.0000 

0.0000 

0.0000 

-4.0000 

0.0000 


_ 0.0000 

0.0000 

12.1943 

0.0000 

0.0000 

0.0000 

0.1890 

-4.0000 

J 


r- 2. 2032 

-2.2032 

-0.8354 

-0.8354 

0.00001 






-0.0840 

-0.0848 

-0.0494 

-0.0494 

0.0000 






-0.0166 

0.01S6 

0.0000 

0.0000 

0.0647 





R„ = 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 






0.0000 

0.0000 

0.0000 

0.0000 

0.0000 






11.2379 

-11.2379 

29.0513 

-29.0513 

9.6847 






-7.8229 

-7.8229 

0.0000 

0.0000 

0.0000 






_ 0.0000 

0.0000 

0.0000 

0.0000 

-6. 6502 J 






*M = 0.8; h 20,000 ft; ct 0 = 1.6°; A = 45° 
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TABLE 2. - GAIN MATRICES* 



V 

a 

8 


e 

P 

q 

r 


r - 0.0519 

- 1.6309 

- 0.0001 

- 0.0432 

2.1180 

- 0.0214 

0.2558 

- 0 . 0239 ] 


- 0.0614 

- 1.0436 

- 0.6015 

0.0630 

1.9146 

0.0475 

0.3378 

0.1859 

lx — 

xp 

- 0.0081 

- 0.4667 

0.1671 

- 0.0406 

0.4541 

- 0.0268 

0.0250 

- 0.0426 

0.0232 

0.6063 

- 0.0928 

0.0333 

- 0.6926 

0.0213 

- 0.0348 

0.0273 


L - 0.0524 

- 0.3968 

- 0.1964 

-0.0009 

0.7165 

- 0.0090 

0.0442 

0 . 1586 . 


V 

a 

8 

♦ 

e 

P 

q 

r 


r 0.0520 

1.4288 

- 0.9317 

0.9512 

02.1145 

0.0205 

- 0.2196 

0 . 0521 ] 


0.0615 

1.2344 

0.1713 

- 0.0608 

- 1.9159 

- 0.0348 

- 0.2036 

- 0.1015 


0.0082 

0.3162 

- 0.0769 

0.0431 

- 0.4550 

0.0225 

- 0.0428 

0.0515 


- 0.0232 

- 0.4681 

0.0552 

- 0.0362 

0.6935 

- 0.0177 

0.0487 

- 0.0424 


l 0.0524 

0.4089 

- 0.0353 

- 0.0013 

- 0.7161 

0.0092 

- 0.0130 

- 0 . 02 56 J 


<5hL 

«hR 

«aL 

*aR 

6R 





f 0.6121 

0.2711 

0.3173 

- 0.5317 

- 0.20121 





0.0015 

0.2122 

- 0.3388 

0.1843 

0.0414 





0.2607 

0.0542 

0.2295 

- 0.2846 

- 0.0907 





- 0.2221 

- 0.0384 

- 0.1941 

0 - 2630 

0.0875 





L - 0.2622 

- 0.1349 

- 0.1708 

0.1445 

0 . 0993 J 





*M = 0.8; h = 20,000 ft; A = 45° 



Fig. 1 Implicit model-following control law 
structure. 



Fig. 2 Explicit model-following control law 
structure. 
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(a) Bank angle. 


Model 

Aircraft 



(b) Xaw rate. 



Model 

Aircraft 


0 1 2 3456789 10 

Time, sec 

(c) Pitch rate. 

Fig. 7 Model -following response to aileron command input. 
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